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FOR MONITORING COASTAL ZONE ENVIRONMENTS^ 

R. A. Weismiller, S. J. Krlstof, D. K. Scholz 
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P. E. Anuta and S. M. Momin 
ABSTRACT 

Development of satisfactory techniques for detecting change in 
coastal zone environments is required before operational monitoring 
procedures can be established. In an effort to meet this need a 
study was directed toward developing and evaluating different types 
of change detection techniques, based upon computer-aided analysis 
of Landsat multispectral scanner (MSS) data, to monitor these 
environments . 

The Matagorda Bay estuarine system along the Texas coast was 
selected as the study area. Land use within this area is divided 
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principally among agriculture, rangeland, urban. Industry, recrea- 
tion and large marsh-covered tracts. 

Four change detection techniques were designed and implemented 
for evaluation: a) post classification comparison change detection, 

b) delta data change detection, c) spectral/temporal change classi- 
fication, and d) layered spectral/ temporal change classification. 
Each of the four techniques was used to analyze a Landsat MSS tem- 
poral data set to detect areas of change of the Matagorda Bay 
region. 

The post classification comparison technique reliably identi- 
fied areas of change and was used as the standard for qualitatively 
evaluating the other three techniques. The layered spectral/tem- 
poral change classification and the delta data change detection 
results generally agreed with the post classification comparison 
technique results; however, many small areas of change were not 
identified. Major discrepancies existed between the post classi- 
fication comparison and spectral/ temporal change detection results. 
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INTRODUCTION 

Our nation's coastal zones are locales of unending change. 
Effective inventorying and monitoring of these changes are required 
for resource management. Conventionally, monitoring of these areas 
has been accomplished either by local on-site surveys or photo- 
interpretation of aerial photographs, ano change detection has re- 
quired the manual comparison of various single date analyses. For 
large area inventories, both local surveys and photointerpretative 
techniques require large data collection efforts and are time con- 
suming and subjective in nature. . 

Consequently, it would be desirable to develop and implement 
a computerized change detection procedure suitable for large area 
inventories. (1) However, before operational monitoring procedures 
can be established, satisfactory computer-aided change detection 
techniques must be developed. In order to meet this need, this 
study was initiated toward developing and evaluating various change 
detection techniques based upon computer-aided analysis of Landsat 
multispectral scanner (MSS) data to monitor coastal zone environ- 
ments . 

STUDY SITE 

A portion of the Matagorda Bay estuarine system located along 
the Texas coast was selected for this study. This area is char- 
acterized by a great diversity in geography, resources, climate, 
waterways and estuaries similar to the entire coastal zone of 


Texas. Also, it is subjected to numerous natural hazards such as 
shoreline erosion, land surface subsidence, stream flooding, hurri- 
cane tidal surges and active surface faulting. Lower elevations of 
the coastal area are covered with marshes and swamps, and the stand- 
plains are chiefly sandy materials. The coastal upland areas con- 
sist mainly of clay and mud materials deposited as sediments during 
the Pleistocene era. Land use within this study site is divided 
principally among agriculture, rangeland, urban, industry, recrea- 
tion and large marsh-covered tracts. (2) 

APPROACH 

Data 

Landsat MSS data collected on November 27, 1972 and February 
25, 1975, were used as the principal data sources for this study. 
These data were geometrically corrected (i.e., rotated, deskewed 
and rescaled) (3) , overlaid and precision registered to ground con- 
trol points selected from appropriate U.S. Geologicnl Survey 
(uses) 7^5 minute topographic quadrangle maps. This procedure pro- 
duced a multidate eight-channel data set at a scale of 1:24,000 
and registered points in the data to their exact ground position. 
Reference data utilized to support the analysis of the Landsat 
data included 1970, 1971 and 1975 color and color infrared aerial 
photography. Geologic Atlases compiled by the Texas Bureau of 
Economic Geology, and Spectral Environmental Classification over- 
lays produced by Lockheed Electronics Corporation, 
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Sln^le Date Analysis 

A non-supervlsed approach utilizing the maximum likelihood 
classifier as implemented by LARSYS (4) was used to produce a clas- 
sification for each date for areas corresponding to individual 
uses 7^5 minute quadrangle maps. The resultant classifications in- 
cluded informational categories o<" agriculture, rangeland, grass- 
lands, wooded areas, swamp timber, fresh and brackish marshes, mud 
and sand flats, beach and beach ridges, residential areas, and 
open water classes. This initial study demonstrate.:! that digital 
analysis of Landsat MSS data could be used to inventory the coastal 
environment effectively. (5) 

Change Detection Analysis 

Four specific test sites were selected within the Matagorda 
Bay estuarine system for Implementation and evaluation of change 
detection techniques. These were the areas included in the Aust- 
well. Port Lavaca-E, Port O'Connor and Pass Cavallo-SW USGS 7^5 
minute topographic quadranglea .. 

Four techniques for detecting change were designed and im- 
plemented for evaluation (6): a) post classification comparison, 

b) delta data change detection, c) spectral/temporal change clas- 
sification, and d) layered spectral/temporal change classification. 
All four techniques require registration of the data prior to 
analysis. A brief description of each technique follows. 



Post Classification Comparison Change Detection (7) 

Post classification change detection is based on the compar- 
ison of independently produced spectral classifications . The com- 
parator is simply a processor that "compares" two standard LARSYS 
classifications utilizing class pairs specified by the analyst and 
generates a results tape indicating areas of change. Thus, the 
change classes derived are defined by the analysts rather than 
being identified spectrally using pattern recognition techniques. 

By properly coding the classification results for times tj^ and t2» 
the analyst can produce change maps which show a complete matrix 
of changes. In addition, selective grouping of classification 
results allows the analyst to observe any subset of changes which 
may be of interest. 

Delta Data Change Detection (8) 

The delta data change detection technique is based upon the 
classification of a multlspectral difference data set. A delta 
(subtraction) transformation combines two n-channel multlspectral 
data sets obtained at different times and produces a multlspectral 
delta data set having n-channels. Since the LARSYS system assumes 
all data samples used are non-negative, a bias is added to each 
difference so that the resulting delta data is non-negative. Thus, 
the complete transformation is simply: 





(tj^) + b 


k 


Where: X = Multlspectral value for channel k; i,j = 1,,.,N 
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assuming an N x N image 
i = Row 
j = Colvrnm 
k = Channel 

b, = Bias for Channel k 
k 

= First date 
t2 - Second date 

This approach assumes that specific changes from one time to 
another will produce a non-negative result which can be detected 
either by examining the images directly or by first classifying 
the delta data, then examining the results displayed in image 
form. The classification method is similar to that for classifi- 
cation of ordinary multispectral data; however, it emphasizes 
classification of multispectral change instead of changes in 
multispectral classification. 

Spectral/Teraporal Change Classification 
This technique is based upon a single analysis of a multidate 
data set to identify change. In a two-date Landsat data set, eight 
channels of data would be analyzed using standard pattern recog- 
nition techniques. By using data sets collected under similar 
conditions at nearly the same day of the year but from different 
years, one would expect that classes where change is occurring 
would have statistics significantly different from where change 


had not occurred and could be so identified. 
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Layered Spectral/Temporal Approach (9,10,11) 

The layered or decision tree classifier is essentially a max- 
imum likelihood classifier using multi-stage decision logic. It 
is characterized by the fact that an unknown sample can be classi- 
fied into a given class using one or more decision functions and 
channels in a successive manner. This classification strategy can 
be most easily illustrated by a tree diagram consisting of a root 
node, a number of non-terminal nodes or decision stages (layers) 
and terminal nodes (i.e,, the decision making procedure terminates 
with the unknown sample being assigned to the class at a terminal 
node). A non-terminal node is an intermediate decision, its imme- 
diate descendant nodes representing the possible outcomes of that 
decision. 

To specify a decision tree uniquely, two sets of information 
are necessary: a) how the terminal and non-terminal nodes are 

linked and b) the decision functions and channels of all the non- 
terminal nodes. The decision tree can be constructed manually or 
by an automatic optimized logic tree design procedure. For the 
purpose of change detection a hybrid of these two procedures was 
utilized. Decision trees for each of the two dates, t^ and t^, 
were obtained automatically and then manually linked introducing 
within the tree a logic for detecting the desired changes. 
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RESULTS 

A quantitative evaluation of the change detection procedures 
developed would require the collection of a comprehensive set of 
reference data (ground observations and/or aerial photography) 
coincident with the Landsat overpass. However, coincident refer- 
ence data were not available for this particular project (i.e., 
reference data for the 1972 Landsat data consisted of aerial pho- 
tography collected in 1970, 1971). Thus, only a qualitative eval- 
uation of the results of the change detection techniques could be 
accomplished. 

Since the comparison of results from two single date classi- 
fications is the basis of the post classification comparison 
change detection procedure, the results of this technique were 
used as the standard for evaluating the results from the other 
three procedures. Correlation of the unsupervised spectral classes 
derived from single date classifications with informational classes 
was accomplished by: 

a) associating the spectral classifications with aerial 
photography when appropriate, 

b) utilizing a ratio, A = ^ calculated for each spectral 

XK. 

class, where V is the relative intensity of the mean 
spectral values of the visible wavelengths {(0.5 to 0.6pm) 
(0.6 to 0.7pm)} and IR is the relative intensity of the 
mean spectral values of the reflective infrared 
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wavelengths {(0.7 to 0.8pm) + (0.8 to 1.1pm)} and 
c) by summing the relative intensity values ("summed re- 
sponse") of all four bands to determine the magnitude 
of relative spectral responses for each spectral class. 

By observing the aerial photography, the ratio A and the summed 
response, the analyst delineated major vegetation and land use 
categories within the coastal zone area. Results from the Port 
O'Connor quadrangle, which are representative of the four test 
areas, are presented for illustrative purposes. 

Utilizing the post classification comparison technique, the 
November 1972 and February 1975 dates for the Port O'Connor test 
site were classified into 25 and 23 spectral classes, respectively. 
These classes were grouped into six ecological categories: 
a) urban, b) woody/herbaceous vegetation, c) submerged vegetation 
and tidal flats, d) spoil areas, e) water and f) burned areas. 

The logic comparator was constructed to identify the class pairs 
listed in Table 1. The map produced by this comparison (Figure 1) 
illustrates the areas of change. 

Cluster analysis of the delta data set produced thirteen 
delta spectral classes representing seven identifiable information- 
al classes (Table 2). These information classes were not in a 1:1 
correspondence with the post classification comparison classes but 
rather were dispersed over a number of the classes. However, 
visual inspection of the delta data change results map showed a 



good general agreement between the post classification comparison 
and delta data change results in that the no-change water/land 
interfaces and the delineation of bum areas were in close agree- 
ment. Pixels classified into the mixed and urban change classes 
were widely dispersed through the land and marsh areas. This 
situation may be related to the inherent difficulty of using 
present Landsat resolution and wavelength bands to discriminate 
reliably among certain surface features and not to the change 
detection procedure. 

To evaluate the spectral/ temporal approach to detecting 
change, three areas from the eight-channel data set were clustered 
to identify potential change classes. These included the agricul- 
tural region, intercoastal waterway and the barrier Island, areas 
where change was expected to be occurring. Statistics developed 
during clustering were used to identify potential change classes 
using the following empirically-derived criteria: 

a change class was probable if the ratio (A = for each 
of the two dates differed by more than 0.30, and the summed 
response differed by more than 25.0. 

Fifty-two spectral classes were identified resulting in five infor- 
mation classes (Table 3), four of which represented change. 

Comparison of the results of the spectral/ temporal approach 
to the post classification comparison technique indicated major 
disagreement. Pir .Is classified into the woody to woody class 


by the post classification comparison technique were classified 
as vegetation to soil, and soil to vegetation change classes 
(Table 3) by the spectral/ temporal change detection technique. 
Similarly, the pixels classified into the sahmerged-'Submerged, any 
non-subraerged-submerged and water-water classes were classified 
as a water to soil change class. Only the burn change cluss 
agreed well with the post classification comparison change results. 

Detection of change using the layered spectral/temporal 
approach Involved developing a layered decision tree based upon 
tlmej^ statistics for time^ (November 1972) and attaching a layer 
decision tree based upon time 2 statistics for time 2 (February 
1975) to the appropriate terminal nodes of tirae^ producing a 
"Change decision tree" . Results of this procedure produced 47 
spectral classes and 18 information classes, 10 of which were 
change classes (Table 4) . 

Results of the layered spectral/ temporal approach were in 
good agreement with those obtained by the post classification 
comparison procedure. Burn areas and the no-change land/water 
interfaces were very similar. Urban areas were well identified; 
however, occasionally the dredged spoil areas along the intra- 
coastal waterway were placed into the urban class due to the 
similarities in spectral response. Some of the confusion between 
classes was attributed to the trimming of the decision trees 
(which was necessary to accomodate the tree in computer memory) . 
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CONCLUSIONS 

In devising a change detection procedure one property to con- 
sider is the relative complexity of the method; all other things 
being equal, a simpler approach would be more desirable for imple- 
mentation reasons. However, more complex methods may have more 
perfomiance potential in the long run. On a scale of increasing 
complexity the four methods investigated may be ordered as follows 

1. Delta data change detection . This method requires only 

a simple subtractior fallowed ’'y a single classification. 

2 . Post classification comparison change detection . This 
method requires two separate classifications followed by 
a logical comparison. 

3. Spectral/ temporal change classification . While this 
method requires only a single classification, it is a 
vastly more complex one, requiring more classes and 
probably more features. 

4 . Layered spectral/ temporal change classification . This 
method involves not only a complex classification but 
also a priori knowledge of the logical interrelationship 
of the classes as well. 

In a change detection analysis, change detection error can 
arise from classification errors at date one, date two or both 
times. To evaluate precisely the results of various change 
detection procedures, comprehensive ground information must be 



- 14 - 


available to Identify any error condition. 

Since adequate reference data for a thorough evaluation was 
not available, the post classification comparison change detection 
results were used as the standard for qualitative evaluation of 
the results from the other three procedures. Areas of major change 
were reliably identified by the post classification comparison 
technique. This is due in part to the fact that the analysis pro- 
cedures required are routine and quite well understood. 

Evaluation of the other three change detection techniques 
indicates the following; 

a) The delta data change detection method may be too simple 
to deal adequately with all the factors Involved in detecting 
change in a natural scene. Too much information may be discarded 
from the data in the subtraction process whereby only the four 
band difference data remains from the two sets of original four 
bands. Images created from the delta data may be quite useful, 
however, in qualitatively assessing change by image interpretation. 

b) The spectral/ temporal and layered spectral/temporal change 
detection methods cannot be ruled out at this point as each appears 
to have undeveloped potential. The layered spectral/temporal 
change detection results in particular showed best agreement with 
the post classification comparison results. More complex methods 
usually require more carefully dratm data inputs together with 
greater user understanding to achieve their potential. 
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Table 1. Class Pairs Defined for the Post Classification 
Comparison Change Detection Technique for the 
Port O'Connor Quadrangle. 


Class No. 

1 

2 

3 

4 

5 

* 6 

* 7 

* 8 
* 9 
*10 
*11 


Informational Classes 
Nov 72 Feb 75 

urban urban 

woody /herbaceous woody /herbaceous 

spoil spoil 

submerged submerged 

water water 

any non-burn burn 

any non-spoil spoil 

any non-woody woody 

any non-water water 

any non-submerged submerged 

any non-urban urban 


* change classes 
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Table 2. 


Class No. 
* 1 
* 2 

* 3 

* 4 

* 5 

6 

7 

8 
9 

10 

11 

*12 

13 


Classes Delineated by Delta Data Change Detection 
Method and Their Corresponding Post Classification 
Comparison Informational Classes for the Port 
O'Connor Quadrangle. 

Informational Classes 
any non-spoil to spoil 
any non-spoil to spoil 

mixed and any non— submerged to submerged 
mixed and any non-urban to urban 
woody to woody and any non-submerged to 
submerged and any non-urban to urban 
woody to woody 

woody to woody and water to water 

woody to woody 

water to water 

water to water 

water to water 

any non-burn to burn 

water to water 


* change classes 


Table 3. Classes Delineated by Spectral/Temporal Change 
Detection Method for Port O'Connor Quadrangle. 


Class No. 
*1 
*2 
*3 
4 

*5 

6 


Informational Classes 
Vegetation to soil 
Soil to vegetation 
Vegetation to burn 
No change land class 
Water soil 

No change water class 


* change class 
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Table 4. Classes Delineated by Layered Spectral/Temporal 
Change Detection Method for the Port O'Connor 
Quadrangle. 


Class No. 

Information 

Classes 

* 1 

woody 

urban 

2 

urban 

urban 

3 

woody 

woody 

* 4 

urban 

woody 

* 5 

woody 

burn 

* 6 

urban 

burn 

7 

water 

water 

* 8 

spoil 

water 

* 9 

submerged 

water 

*10 

water 

spoil 

11 

spoil 

spoil 

*12 

submerged 

spoil 

*13 

water 

submerged 

*14 

spoil 

submerged 

15 

submerged 

submerged 

16 

water 

confusion 

17 

spoil 

confusion 

18 

submerged 

confusion 


* change classes 


